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Summary 

The suitability of the MARC code for the analysis of 
rotating composite blades was evaluated using a four-task 
process. A nonlinear displacement analysis and subsequent 
eigenvalue analysis were performed on a rotating spring-mass 
system to ensure that displacement-dependent centrifugal 
forces were accounted for in the eigenvalue analysis. Normal 
modes analyses were conducted on isotropic plates with 
various degrees of twist to evaluate MARC’s ability to handle 
blade twist. Normal modes analyses were conducted on flat 
composite plates to validate the newly developed coupled 
COBSTRAN-MARC methodology. Finally, normal modes 
analyses were conducted on four composite propfan blades 
that were designed, analyzed, and fabricated at NASA Lewis 
Research Center. Results were compared with experimental 
data. The research documented herein presents MARC as a 
viable tool for the analysis of rotating composite blades. 


Introduction 

Since the advent of the Advanced Turboprop Project (ATP), 
NASA Lewis Research Center has become the lead center in 
the design and analysis of rotating composite propfan blades 
(ref. 1 ). During the ATP program, the Lewis research commu- 
nity expended considerable effort developing the methodolo- 
gies currently used to structurally analyze composite propfan 
blades. These methodologies can be, or have been, applied to 
other rotating components, such as compressor blades and 
turbine blades. As a result of this effort, the COBSTRAN 
(Composite Blade STRuctural ANalyzer) preprocessor code 
was developed for the generation of composite blade finite 
element models (ref 2). Given blade geometry and material 
selection, the COBSTRAN preprocessor combines composite 
micromechanics and classical laminate theory with a data 
base of fiber and matrix properties to generate a finite element 
model with anisotropic material properties. The COBSTRAN 
code, coupled with MSC or COSMIC NASTRAN, has been 
the primary tool used thus far for the analysis of these blades. 

Recently, the COBSTRAN code was combined with the 
MARC finite element program. The MARC finite element 
package is formulated principally for nonlinear analyses 


(ref. 3). By providing an alternative to NASTRAN, the cou- 
pling of MARC with COBSTRAN has the potential to provide 
the analyst with greater versatility. 

This report documents the research done to evaluate the 
suitability of MARC in the analysis of rotating composite 
blades. The evaluation process consisted of four tasks. First, 
a nonlinear displacement analysis and subsequent eigenvalue 
analysis were performed on a simple spring-mass system 
subjected to a centrifugal load. This was done to ensure that 
displacement-dependent centrifugal forces are included in the 
stiffness matrix used during an eigenvalue analysis. Second, 
normal modes analyses were conducted on isotropic plates 
with various degrees of twist, and the results compared with 
experimental data, to evaluate MARC’s ability to handle 
blade twist. Third, normal modes analyses were conducted on 
flat composite plates, and the results compared with experi- 
mental data. The coupled COBSTRAN-MARC methodology 
that was developed in this effort was thus validated. Finally, 
as a direct application of the methodology, normal modes 
analyses were conducted on four composite propfan blades 
that were designed, analyzed, and fabricated at NASA Lewis 
Research Center. Results were compared with holographic 
data obtained from bench tests conducted on these blades. 


Symbols 

a plate length 
b plate width 
D plate flexural rigidity 
E modulus of elasticity 
G shear modulus 
h plate thickness 
/ mass moment of inertia 
K spring stiffness 
L Lagrangian function 
m mass 

r initial spring length 
V potential energy 
y shear strain 


l 


I 


£ direct strain 
k kinetic energy 
A frequency parameter 
jj radial displacement 

v Poisson’s ratio 

p mass density 

<7 direct stress 

T shear stress 

0 twist angle (root to tip) 

Q rotational velocity 
co frequency 


Subscripts: 

x longitudinal 
y transverse 
z normal 


Centrifugal Softening in Eigenvalue 
Analysis 

When conducting a modal analysis on a rotating, flexible 
structure, two types of forces are considered: forces due to the 
component of the centrifugal load that is displacement depen- 
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Figure I. — Sample MARC input deck 


dent and forces due to gyroscopic (or Coriolis) effects. These 
forces are not always accounted for in the respective solution 
sequences of available finite element packages. With regard to 
the dynamic analysis of rotating propfan blades, gyroscopic 
effects are considered negligible because the angular velocity 
is generally lower than the first modal frequency of the 
structure (refs. 4 and 5). However, forces due to the compo- 
nent of centrifugal force that is displacement dependent are 
very important (ref. 6). 

The initial step in evaluating MARC’s suitability for ana- 
lyzing rotating, flexible blades was to execute a nonlinear 
displacement analysis and a subsequent eigenvalue analysis 
on a rotating spring-mass system. This simple test was used to 
determine if the component of centrifugal force, dependent 
upon the translational degrees of freedom lying in the plane of 
rotation, was accounted for in the stiffness matrix during the 
eigenvalue analysis. The terms in the stiffness matrix that 
account for these displacement-dependent centrifugal loads 
are referred to as “centrifugal-softening” terms. 

Centrifugal loads are displacement dependent. As a result, 
to dynamically analyze a finite element model of a flexible 
structure experiencing centrifugal loading, a nonlinear dis- 
placement analysis is initially required. The nonlinear prob- 
lem is solved by a series of linear analyses which update load 
and stiffness matrices. Subsequent to the nonlinear displace- 
ment analysis, a modal analysis is performed with the updated 
mass and stiffness matrices that result from the nonlinear 
displacement analysis. 

By modeling a simple spring-mass system in the plane of 
rotation, MARC displacement and frequency results at speed 
were easily compared with the displacement and frequency 
results obtained from the equation of motion derived in appen- 
dix A. For example, when the followingnondimensional values 
are assumed: spring stiffness K of 500, initial spring length 
v of 1, end mass m of 1 , and rotational velocity Q of 10 rad/sec, 
the differential equation of appendix A will yield a steady state 
displacement of 0.25 and a frequency of 20 rad/sec. MARC’s 
nonlinear analysis should converge upon these values. If the 
stiffness matrix does not account for the centrifugal-softening 
terms, a frequency of 22.36 rad/sec will result. 

Modeling a simple spring-mass system with MARC is not 
straightforward. There is no provision for defining a point 
mass at the end of a spring. To resolve this problem, the end 
mass was modeled with a MARC finite element of extremely 
small proportions such that the centroid of the element was 
infinitesimally close to the end node of the spring. The model 
then consisted of a spring defined between two nodes: one 
node constrained in all degrees of freedom, and the other 
coupled to the small finite element representing the end mass. 
The finite element geometry and material density were defined 
so that the mass of the element was equal to 1. 

Figure 1 shows the MARC input deck required to run 
nonlinear displacement analysis and subsequent eigenvalue 
analysis. The MARC input deck is divided into three sections: 
the parameter definition section, the model definition section, 


and the history definition section (ref. 7). In the parameter 
definition section, the DYNAMIC card flags the eigenvalue 
analysis, whereas the nonlinear displacement analysis is flagged 
by the FOLLOW FOR and LARGE DISP cards (ref. 8). The 
model definition section includes all geometry and material 
property data, boundary conditions, and initial loading condi- 
tions. In particular, the centrifugal force is applied to the 
model by specifying an IBODY load type of 100 (ref. 9), along 
with the rotational velocity squared, within the DIST LOAD 
card block (ref. 8). The rotational axis is specified with the 
ROTATION A card block (ref. 8). The history definition 
section controls the loading throughout the subsequent incre- 
ments of the analysis. The DIST LOAD card blocks included 
in this section are used only to update the stiffness matrix 
before an eigenvalue analysis is conducted; therefore, rota- 
tional velocity squared is not specified on these cards. 

The element library of the MARC code contains over 100 
different element formulations. Only shell elements were 
chosen for evaluation. Among the available shell elements, 7 
elements looked practical for rotating blade applications: 
elements 3, 22, 26, 49, 50, 72, and 75 (ref. 9). The end mass 
of the spring-mass system was modeled with each of these 
elements in turn. Analyses conducted with elements 3, 22, 26, 
and 75 were found to correctly solve for the displacement and 
eigenvalue of the spring-mass problem previously discussed. 
Appendix B summarizes the formulations for these elements. 
Table I gives the results of analyses using these four elements. 
It is not understood why the remainder of the shell elements 
tested did not correctly solve the spring-mass problem. 

When modeling a rotating, flexible structure with element 
3, 22, 26 or 75, MARC will correctly perform a nonlinear 
displacement analysis and subsequent eigenvalue analysis 
without special programming to update the stiffness matrix 
with centrifugal-softening terms. Because of the two- 
dimensional nature of elements 3 and 26, the remainder of this 
study was restricted to elements 22 and 75. 


TABLE I. -RESULTS FROM ROTATING SPRING-MASS 
SYSTEM PROBLEM 


MARC element 

Displacement 

Converged in 
increment 

Eigenvalue 

3 

□ 

0.25129 

1 

19.9996 

3 

A 

.25089 

1 

19.9995 

22 

□ 

.25129 

2 

19.9996 

22 

A 

.25089 

2 

19.9995 

26 

□ 

.25129 

2 

19.9996 

26 

A 

.25089 

2 

19.9995 

75 

□ 

.25129 

1 

19.9996 

75 

A 

.25089 

1 

19.9995 
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To demonstrate MARC’s ability to properly handle the 
physics of centrifugal loading during a modal analysis run, a 
flat plate was modeled and analyzed at rotational velocities of 
2000, 4000, and 5000 rpm for plate setting angles of 0° and 
90° with respect to the plane of rotation. When the plate is 
lying in the plane of rotation (0° orientation), the stiffness 
associated with the translational degrees of freedom perpen- 
dicular to the plate does not experience centrifugal softening; 
however, when the plate is perpendicular to the plane of 
rotation (90° orientation), the stiffness associated with the 
same degrees of freedom does experience centrifugal soften- 
ing. Therefore, it was expected that the first bending fre- 
quency of the flat plate lying in the plane of rotation would be 
higher at rpm than the first bending frequency of the flat plate 
perpendicular to the plane of rotation. 

The flat plate modeled was 7 in. long, 3 in. wide, and 1/8 in. 
thick. The finite element model of the flat plate was composed 
of 65 quadrilateral elements (because of its simplicity, MARC 
element type 75 was used). Figure 2 displays the model and 
lists the material properties used. All six nodal degrees of 
freedom were constrained along the short edge of the plate 
closest to the axis of rotation. The axis of rotation was 7 in. 
from the constrained edge of the plate. 

Figure 3 shows the first bending frequencies of the flat plate 
at rpm for the two plate orientations of 0° and 90°. As 
expected, at rpm the first bending frequency of the flat plate 
lying in the plane of rotation was higher than that for the flat 
plate perpendicular to the plane of rotation. Identical results 
were obtained by analyzing an equivalent finite element 
model subjected to the same boundary conditions with MSC/ 
NASTRAN using extensive DMAP programming (ref. 6). 
With MARC, however, no special programming was required 
to add centrifugal-softening terms into the global stiffness 
matrix utilized by the modal analysis. 

Note that several of the MARC element types (elements 50 
and 72) that failed the rotating spring-mass problem were 
again tested in this flat plate analysis. The frequencies were 
identical to those shown in figure 3. Erroneous results with 
centrifugal loading problems occurred when the finite ele- 
ment model was constrained with springs. Because situations 
may arise when spring elements are used to simulate the hub 
stiffness associated with a particular rotating blade (ref. 10), 
it is recommended that finite element models of rotating 
structures be created only with the plate elements that prop- 
erly solved the rotating spring-mass problem discussed here. 


Isotropic Twisted Plate Study 

Blade twist is one geometric aspect of rotating blades that 
concerns the analyst. When modeling twisted blades with 
quadrilateral elements, element nodes generally are notcopla- 
nar. For certain finite elements, this may result in an additional 
stiffness that is inaccurate (ref. 11). Reference 12 presents 



3 in. 


Properties 

Modulus of elasticity = 1 5.5x1 0 6 psi 
Shear modulus = 6.5x1 0 6 psi 
Density = 0.163 lb/in. 3 


1/8-in. thickness 


Figure 2. — Finite clement model of flat plate. 



Rotational velocity, rpm 

Figure 3. — Effect of centrifugal softening on first bending frequencies of 
flat plate. 


experimental frequencies for various cantilevered isotropic 
twisted plates of uniform thickness. A normal modes analysis 
was conducted on a selected number of twisted plates and 
compared with the experimental results found in reference 1 2. 
This was done to assess MARC’s ability to perform such 
analyses when blades are modeled with element 75 or 22. 
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Reference 12 presents experimental frequencies and mode 
shapes for various cantilevered twisted plates machined out of 
7075-T6 cold-drawn aluminum bar stock (fig. 4). Refer- 
ence 13 provides detailed information on the experimental 
methods used. Twenty twisted plates, consisting of all com- 
binations of two aspect ratios ( alb - 1 and 3), two thickness 
ratios (b/h = 5 and 20), and five twist angles (0 = 0°, 15°, 30°, 
45°, and 60°) were machined. All the plates had a rectangular 
cross section with a width b of 2 in. Frequencies documented 
in reference 1 2 were reported in the form of the standard plate 
nondimensional frequency parameter: 


i 2 j P X 

A = coxa x J- — 

V D 

where D is the plate flexural rigidity, 


12x(l- v 2 ) 

This study was restricted to plates with thickness ratios of 
20 and aspect ratios of 3 because these plates reflect the 
geometric traits of most propfan blades. Two sets of finite 
element models were generated for plates with twist angles of 
0°, 30°, and 60° (fig. 5): one set comprising 384 4-noded 
quadrilateral elements (MARC element 75), and one set 
comprising 60 8-noded quadrilateral elements (MARC 
element 22). For each model, nodal degrees of freedom were 
fully constrained along one of the short edges. A modulus of 
elasticity E of 10.3X10 6 psi, Poisson’s ratio v of 0.33, and 
mass density p of 0.101 lb/in. 2 were assumed. 

A MARC normal modes analysis was conducted on each 
of the modelsgenerated. Analytical mode shapes and frequen- 
cies were recorded for the first three bending modes, the first 
two torsional modes, and the first edgewise mode. All fre- 
quencies were put in the form of the standard plate 


2 




Figure 4. — Twisted cantilevered plate. 


nondimensional frequency parameter X. Results from the 
analyses are presented in figures 6 and 7 for the models 
consisting of elements 75 and 22, respectively. Table II pre- 
sents a comparison between analytical and experimental 
nondimensional frequencies. An experimental first edgewise 
frequency is not available for the plate with 0° twist. Gener- 
ally, analytical frequencies correlated very well with experi- 
mental frequencies. However, poor correlation was noted for 
the 3rd bending mode associated with the plate with 30° twist, 
and for all edgewise modes. 

Flat Composite Plate Analysis 

Composite materials are being used more frequently in the 
construction of rotating blades. However, few finite element 
codes provide composite analysis capability that reflects 
actual composite blade construction. As mentioned in the 
introduction, the COBSTRAN preprocessor code was devel- 
oped for the generation of composite blade finite element 
models (refs. 2 and 14). Formerly, COBSTRAN generated 
either MSC or COSMIC NASTRAN input decks. The 
COBSTRAN code has recently been modified so that a 
MARC input deck for a blade model can also be created. At 
this time, only a model using MARC element 75 is output by 
COBSTRAN. In order to verify the modified version of 
COBSTRAN, as well as MARC’s ability to analyze compos- 
ite blades, a normal modes analysis was performed on flat 
composite plates that were previously dynamically bench 
tested. MARC eigenvalues were compared with those gener- 
ated by MSC/NASTRAN and the dynamic bench test. Three 
flat graphite-epoxy composite plates, designated form 1, 
form 2, and form 3, were modeled and analyzed. Figure 8 
shows the finite element model, material properties, geom- 
etry, and ply lay-up for the three plates. Each plate was made 
of 18 plies resulting in a thickness of 0.09 in. The composite 
plies were symmetrical about the midplane of the plates and 
were stacked from the surface to the point of midthickness. All 
plies were oriented with respect to the spanwise axis (from 
base to tip). Starting with the outside ply and working toward 
the point of midthickness, the plates were fabricated with the 
ply groups shown in figure 8. Finite element models of these 
flat plates were generated with COBSTRAN and analyzed 
with both MARC and MSC/NASTRAN. Quadrilateral ele- 
ments were used because, at this time, the COBSTRAN code 
does not provide for the output of a MARC model that has 
triangular elements. The model was constrained in all degrees 
of freedom at the eight nodes along the base of the plate. 

The experimental data used for comparison with the ana- 
lytical results were obtained from a pending reference publi- 
cation (T.J.Sutliffetal., to be published). In this investigation, 
vibration mode frequency and shape evaluations were made 
on three cantilevered composite plate specimens described in 
figure 8. 
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0° Twist 



Figure 5. — Twisted plate study; MARC finite clement models. 
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Figure 6. — Twisted plate study; MARC clement 75 results. (Nondimensional frequencies are 
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Figure 7. — Twisted plate study; MARC element 22 results. (Nondimensional frequencies arc 
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presented.) 
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TABLE II. -TWISTED PLATE STUDY 


Mode 

Experiment 

MARC 
element 75 
QUAD-4 

Difference 

from 

experiment, % 

MARC 
element 22 
QUAD-8 

Difference 

from 

experiment, % 

Twist angle, 0° 

1st Bending 

3.32 

3.40 

2.41 

3.14 

-5.42 

2nd Bending 

20.76 

21.30 

2.60 

19.70 

-5.11 

1st Torsion 

19.58 

20.61 

5.26 

19.93 

1.79 

3rd Bending 

58.65 

60.20 

2.64 

55.71 

-5.01 

2nd Torsion 

62.43 

64.90 

3.96 

60.08 

-3.76 

1st Edgewise 

N/A 

61.70 

N/A 

58.34 

N/A 


Average 3.37 

Average —3.50 

Twist angle, 30° 

1st Bending 

3.28 

3.40 

3.66 

3.14 

-4.27 

2nd Bending 

18.03 

18.94 

5.05 

17.45 

-3.22 

1st Torsion 

23.53 

24.90 

5.82 

24.08 

2.34 

3rd Bending 

46.29 

56.01 

21.00 

52.02 

12.38 

2nd Torsion 

71.61 

77.25 

7.88 

72.97 

1.90 

1st Edgewise 

57.98 

71.22 

22.84 

67.50 

16.42 


Average 1 1 .04 

Average 4.26 

Twist angle, 60° 

1st Bending 

3.22 

3.40 

5.59 

3.15 

-2.17 

2nd Bending 

14.09 

15.03 

6.67 

13.74 

-2.48 

1st Torsion 

32.29 

34.18 

5.85 

33.12 

2.57 

3rd Bending 

47.68 

5 1 .02 

7.01 

47.22 

-0.96 

2nd Torsion 

94.02 

99.80 

6.15 

96.43 

2.56 

1st Edgewise 

73.47 

85.87 

16.88 

81.88 

1 1 .45 


Average 8.02 

Average 1 .83 



8 in. 


Material 

Fiberite 21 48A1 A epoxy 
Properties 

Fiber volume ratio 0.604 

Weight density, lb/in. 3 0.560x1 0 -1 

Longitudinal modulus, lb/in. 2 0.347x1 0 8 

Transverse modulus, lb/in. 2 0.1 1 0x1 0 7 

Shear modulus, lb/in. 2 0.800x1 0 6 

Poisson's ratio 0.300 

Composite ply strengths 

Longitudinal tensile, lb/in. 2 0.1 46x1 0 6 

Longitudinal compressive, lb/in. 2 .... 0.1 46x1 0 6 

Transverse tensile, lb/in. 2 0.1 50x1 0 5 

Transverse compressive, lb/in. 2 0.21 0x1 0 5 

Intralaminar shear, lb/in. 2 0.1 1 8x1 0 2 


Ply lay-up 
Form 1 (0 
Form 2 (0 
Form 3 (0 

t 

Plate 

surface 


0 90 90 

0 45 90 

45 0 -45 


0 0 90 

-45 0 45 

0 45 0 


90 

0) 

0 

-45) 

-45 

0) 


t 


Plate 

midthickness 


0.09-in. thickness 


Figure 8. — Composite flat plate model, forms 1, 2, and 3. 
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2918.0 Hz 


1 1 


I 


Figure 9. — Results for composite plate, form 1 . 
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Figure 10.— Results for composite plate, form 2. 
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Not Not 

available available 



90.1 Hz 


496.9 Hz 


561.3 Hz 1555.5 Hz 


1590.6 Hz 




Not 

available 


88.4 Hz 

523.4 Hz 

550.4 Hz 


1644.4 Hz 

1545.1 Hz 


89.2 Hz 

466.8 Hz 

535.7 Hz 

1455.6 Hz 

1465.9 Hz 


1493.4 Hz 


Figure 11.— Results for composite plate, form 3. 
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TABLE III. -SUMMARY OF COMPOSITE FLAT PLATE ANALYSES 


Mode 

Experiment 

MARC 

Difference 

from 

experiment, % 

NASTRAN 

Difference 

from 

experiment, % 

Form 1 

1st Bending 

91.6 

83.9 

-8.4 

91.9 

0.3 

2nd Bending 

554.3 

524.8 

-5.3 

572.8 

3.3 

3rd Bending 

1526.3 

1465.2 

-4.0 

1595.5 

4.5 

4th Bending 

2918.0 

2860.8 

-2.0 

3106.3 

6.4 

Form 2 

1st Bending 

89.0 

81.8 

-8.1 

90.1 

1.2 

1st Torsion 

415.7 

518.7 

24.8 

441.8 

6.3 

2nd Bending 

539.1 

511.4 

-5.1 

564.3 

4.7 

2nd Torsion 

1318.8 

1621.6 

23.0 

1409.0 

6.8 

3rd Bending 

1491.8 

1432.8 

-4.0 

1585.1 

6.2 

Form 3 

1st Bending 

89.2 

88.4 

-0.9 

90.1 

1.0 

1st Torsion 

466.8 

523.4 

12.1 

496.9 

6.4 

2nd Bending 

535.7 

550.4 

2.7 

561.3 

4.8 

Mixl a 

1455.6 

N/A 

N/A 

1555.5 

6.9 

2nd Torsion 

1465.9 

1644.4 

12.2 

N/A 

N/A 

3rd Bending 

N/A 

1545.1 

N/A 

N/A 

N/A 

Mix2 b 

1493.4 

N/A 

N/A 

1590.6 

6.5 


"First coupled bending/torsion. 
h Sccond coupled bending/torsion. 


Figures 9 to 1 1 display mode shapes and eigenvalues from 
the MARC and NASTRAN analyses, along with the corre- 
sponding test results, for each of the three forms of the 
composite plate. Table III presents the tabulated results. From 
this study, some general observations can be made. First, 
MARC torsional frequencies are generally higher than 
NASTRAN or experimental torsional frequencies, whereas 
MARC bending frequencies are generally lower than 
NASTRAN or experimental bending frequencies. The second 
bending mode of form 3 is an exception. Second, MARC 
eigenvalues correlate better with experimental results when 
the ply pattern is recurrent throughout the laminate. For 
example, the composite ply lay-ups for forms 1 and 3 are 
recurrent, whereas the composite ply lay-up for form 2 is 


irregular. Third, MARC-generated mode shapes tend to be 
very symmetrical, whereas NASTRAN -generated mode shapes 
tend to be unsymmetrical. With regard to the composite plate 
form 3 analysis (fig. 11), note that certain mode shapes were 
missed by MARC, and others, by NASTRAN. 

Analysis of Composite Propfan Blades 

The NASA Lewis Research Center designed, analyzed, and 
fabricated counter-rotating composite propfan blades for a 
0.55-scale cruise missile wind tunnel model (ref. 15). Fig- 
ure 12 shows a drawing of this cruise missile. Two sets of 
forward and aft blades were designed and fabricated: one set 
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Figure 12. — Cruise missile conceptual design. 


C-90-1 4578 



Figure 13. — CM 1 propfan blade. 


is designated the CM1 series (fig. 13), and the other set is 
designated the CM2 series (fig. 14). Details on the geometric 
design of these blades can be found in references 16 and 17. 
Originally, these blades were dynamically analyzed with 
MSC/NASTRAN's version 65C (refs. 18 and 19). Recently, a 
normal modes analysis was conducted on these propfan blades 
with MARC. This study was carried out as the final step in the 
evaluation of MARC for use in analyzing rotating composite 
propfan blades. Mode shapes and frequencies resulting from 
the analyses were compared with those obtained from MSC/ 
NASTRAN and with those obtained from holographic testing 
(ref. 20). 

All four propfan blades were fabricated with graphite- 
epoxy composite plies. Each ply had a thickness of 0.0032 in. 
Table IV presents the material properties for these plies. The 


TABLE IV - MATERIAL PROPERTIES FOR 
GRAPHITE EPOXY PLIES 


Fiber volume ratio 0.604 

Weight density, Ib/in. 3 0.560.x I0~ 1 

Longitudinal modulus, Ib/in. 2 0. 194x10 s 

Transverse modulus, Ib/in. 2 0.120X10 

Shear modulus, Ib/in." 0.700 X 10*’ 

Poisson's ratio 0.310 

Longitudinal tensile strength, Ib/in. 2 0.266 X ID*’ 

Longitudinal compressive strength, Ib/in." 0.266 X I0 6 

Transverse tensile strength, Ib/in. “ 0.930X10 4 

Transverse compressive strength, Ib/in." 0.930 X It) 4 

Intralaminar shear strength, Ib/in. - 0.130X10" 
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Figure 14. — CM2 propfan blade. 
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composite plies were symmetrical about the midplane of the 
blades and were stacked from the surface to the point of mid- 
thickness. Because the thickness of the blades varies from root 
to tip and from leading edge to trailing edge, the composite 
plies were contoured such that ply size diminished when 
working from the surface to the point of midthickness. All 
plies were oriented with respect to the spanwise/stacking axis 
(from base to tip; see fig. 15). Angles were positive when plies 
were oriented toward the trailing edge of the blades. Starting 
with the outside ply and working toward the point of 
midthickness, the CM1 blades were fabricated with a repeat- 
ing ply group of [0°, 0°, 45°, 0°, 0°, -45°]. The CM2 blades 
were fabricated with a repeating ply group of [20°, 20°, 65°, 
20°, 20°, -25°]. Reference 21 provides further details on the 
fabrication of these propfan blades. 

All MARC finite element models pertaining to the cruise 
missle propfan blades were generated with the COBSTRAN 
preprocessor code and utilized MARC’s element 75. When 
generating MARC finite element models of composite propfan 
blades, equivalent material properties of the composite lami- 
nate are specified for each finite element with MARC’s 
ORTHOTROPIC material definition block (ref. 8). The aniso- 
tropic material properties were defined in a coordinate system 
whose primary axis lies along the stacking axis shown in 
figure 15. The material properties were then incorporated into 
stress-strain laws for the generation of a “compliance” matrix 
(ref. 7) for elements 22 and 75: 
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MARC’scompliance matrix is not fully populated because the 
orthotropic material properties are defined along the material 
axis of symmetry. Thus, the longitudinal and transverse strains 
are uncoupled from the shear stresses, and the shear strains are 
uncoupled from the longitudinal and oansverse stresses. 

As mentioned earlier, the composite ply lay-up for the CM2 
forward and aft blades was a repeating ply group of [20°, 20°, 


65°, 20°, 20°, - 25°] with respect to the stacking axis. Because 
of the simplistic nature of the ply lay-up, it can be seen that 
the axis of material symmetry was off by 20° from the stacking 
axis of the blade. Because COBSTRAN always generates 
finite element models such that the primary axis of the 
material coordinate system lies along the “stacking” axis of 
the blade, and MARC requires properties about the material 
axis of symmetry, COBSTRAN would not directly generate a 
MARC finite element model that accurately reflected the 
CM2 design. This was not a problem with the CM1 design 
because the axis of material symmetry coincided with the 
stacking axis. 

To remedy the problem, the equivalent laminate material 
properties for each element of the CM2 blade were redefined 
within the MARC input deck: specifically, the properties were 
written in a local material coordinate system where the pri- 
mary axis coincided with the material axis of symmetry. This 
was accomplished by using COBSTRAN to regenerate the 
finite element model with a repeating ply group of [0°, 0°, 45°, 
0°, 0°, — 45°] with respect to the stacking axis. By adding 20° 
to the orientation angles defined in the ORIENTATION defi- 
nition data block (ref. 8), the material coordinate system for 
each element was redefined such that the primary axis was 
rotated 20° from the stacking axis about the local element 
normal. The COBSTRAN-generated MARC model then 
accurately reflected the CM2 blade design. 

A normal modes analysis was conducted on all four propfan 
blade models using MARC. MARC finite element models of 
the CM1 forward and aft blades comprised 247 quadrilateral 
elements. In all, the elements were defined by 280 coordi- 
nates: 14 coordinates for each of 20 spanwise stations of the 
blade. CM2 forward and aft models comprised 228 quadrilat- 
eral elements. The elements in the CM2 models were defined 
by 260 nodes: 1 3 coordinates for each of 20 spanwise stations. 
Translational and rotational degrees of freedom of nodes 
closely approximating the location of the blade-to-hub inter- 
face were fully constrained: four nodes at the base of the CM 1 
models, located about the midchord, were fully constrained; 
while three nodes at the base of the CM2 models were fully 
constrained. For comparison, the finite element models docu- 
mented in references 18 and 19 were reinvestigated with 
MSC/NASTRAN’s normal modes analysis using the same 
boundary conditions imposed upon the MARC models shown 
in figure 15. The only difference between the MARC and 
NASTRAN models was the type of elements used: triangular 
elements were used in the NASTRAN models, whereas quad- 
rilateral elements were used in the MARC models. 

Figures 16 to 19 show the resulting mode shapes and 
frequencies obtained from these analyses. MARC and 
NASTRAN mode shapes and frequencies correlated well 
with experimental results. With the exception of the first 
bending frequencies for the CM2 forward and aft blades, 
MARC analytical frequencies were lower than NASTRAN 
analytical frequencies. 
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Figure 1 6.— Results for CM 1 forward blade. 


2881.5 Hz 


IS 






Experimental 


NASTRAN 


MARC 








Figure 1 7.— Results for CM 1 aft blade. 
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Figure 18.— Results for CM2 forward blade. 
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Figure 1 9. — Results for CM2 aft blade. 
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Conclusions 

For the shell elements used in this study, displacement- 
dependent centrifugal forces are accounted for by MARC in 
the dynamic analysis of rotating, flexible structures. Unlike 
NASTRAN (ref. 6), MARC does not require special program- 
ming to account for these forces. However, when springs are 
coupled to a MARC finite element model subjected to cen- 
trifugal loading, a simple test, such as the rotating spring-mass 
problem, should be made on the chosen finite element to 
ensure inclusion of centrifugal-softening terms in the stiffness 
matrix. 

The isotropic twisted plate study and the analyses con- 
ducted on the composite propfan blades show that MARC’s 
quadrilateral elements 75 and 22 are capable of handling the 
complex geometric attributes of blades (i.e., twist and curva- 
ture). The use of equivalent orthotropic material properties of 
composite laminates is suitable when generating MARC mod- 
els of composite blades with COBSTRAN. However, these 
material properties must be defined about the material axis of 
symmetry. A recurrent composite ply pattern is recommended 
for use in the design of the propfan blade. 

Keeping in mind the aforementioned limitations, the use of 
MARC coupled with COBSTRAN offers the analyst an alter- 
nate approach for analyzing rotating composite blades. The 
COBSTRAN-M ARC methodology can be used in place of, or 
in conjunction with, the current COBSTRAN-NASTRAN 
approach. This methodology provides the analyst with greater 
versatility because of MARC’sextensive nonlinear capabilities. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 21, 1992 
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Appendix A 

Equation of Motion for a Rotating Spring-Mass System 



m Mass 

K Spring stiffness 
ft Rotational velocity 

r Initial spring length 

/i Radial displacement 


The Lagrangian L is then 


L = k— V 

= 1/2 m(r+n) 2 £2 2 + 1/2 m f-l/lKn 2 


Applying Lagrange’s equation 


The moment of inertia / about the rotational axis is written as 
l=m(r + n) 2 

The kinetic energy k is written as 

k=\/21Q 2 + 1/2w/} 2 
= M2m(r + n) 2 Q 2 + 1/2 m /j 2 
The potential energy V is written as 
V= 1/2AT/J 2 


dL d_(dL ' 

d/J dt{d/j y 


the equation of motion is derived as 

m(r + /j)C2 2 - Kju -mj} = 0 


m'\i + Kfj = mQ(r + ju) 
m ft + [K - m Q 2 = mQ 2 r 

V ”T" 

centrifugal -softening term 



Appendix B 

Summary of MARC Shell Element Formulations 


Element 3 

Quadrilateral or triangular 

Element 26 

Quadrilateral or triangular 


4-node 


8-node 


2 degrees of freedom 


2 degrees of freedom 


Isoparametric 


Isoparametric 


Plane stress 


Plane stress 


Uniform thickness at nodes 


Uniform thickness at nodes 

Element 22 

Quadrilateral or triangular 

Element 75 

Quadrilateral or triangular 


8-node 


4-node 


6 degrees of freedom 


6 degrees of freedom 


Transverse shear effects 


Transverse shear effects 


Variable thickness at nodes 


Variable thickness at nodes 
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